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We present and discuss our recent study of an eikonal two channel model, in which we reproduce 
the soft total, integrated elastic and diffractive cross sections, and the corresponding forward dif- 
ferential slopes in the ISR-Tevatron energy range. Our study is extended to provide predictions at 
the LHC and Cosmic Rays energies. These are utilized to assess the role of unitarity at ultra high 
energies, as well as predict the implied survival probability of exclusive diffractive central production 
of a light Higgs. Our approach is critically examined so as to estimate the margins of error of the 
calculated survival probability for diffractive Higgs production. 
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I. INTRODUCTION 

The search for unambiguous s-channel unitarity signa- 
tures in ultra high energies soft hadronic scattering, is 
two folded: 

On the one hand, this is a fundamental issue on which 
we have only limited information from the ISR-Tevatron 
experiments. The only direct indication we have on the 
importance of unitarity considerations, derives from the 
observation that soft diffraction cross sections, essentially 
SD (single diffraction), have a much milder energy depen- 
dence than the seemingly similar, elastic cross sections. 
Enforcing unitarity constraints is a model dependent pro- 
cedure. Thus, reliable modeling is essential for the exe- 
cution of our study, leading to predictions of interest for 
LHC and AUGER experiments. 

On the other hand, unitarity considerations in soft scat- 
tering are instrumental for the assessment of inelastic 
hard diffraction rates, specifically, diffractive Higgs pro- 
duction at the LHC. Preliminary information on the im- 
portance and method of this calculation has been ac- 
quired in the study of hard diffractive di-jets at the 
Tevatron[3|, leading to first generation estimates of the 
corresponding survival probabilities. 

This presentation is based on our recent paper 
which utilizes the GLM model [J, |3, 0, d, [1, 0] where 
we numerically solve the s-channel unitarity equation 
in an eikonal model. Our updated results, in the ISR- 
Tevatron range, were obtained from an improved two 
channel model calculations. The specific objectives of 
our study, based on the above, were: 

1) To reproduce the total, integrated elastic and diffrac- 
tive cross sections and corresponding forward differen- 
tial slopes in the ISR-Tevatron energy range, and to ob- 
tain predictions for these observables at LHC and Cosmic 
Rays energies. 

2) To calculate the survival probabilities of inelastic hard 
diffractive processes [1, Q. This requires precise knowl- 
edge of the soft elastic and diffractive scattering ampli- 
tudes of the initial hadronic projectiles. As we noted, 
it is of particular importance for the assessment of the 



discovery potential for LHC Higgs production in an ex- 
clusive central diffractive process. 

3) Some of the fundamental consequences of s-channel 
unitarity in the high energy limit are not clear, as yet. 
We examin the approach of the scattering amplitudes to 
the black disc bound. 

4) We estimate the margin of error of our predicted sur- 
vival probabilities, based on a critical analysis of our 
model. 



II. THE GLM MODEL 

The main assumption of the two channel GLM model 
is that hadrons are the correct degrees of freedom at 
high energies, diagonalizing the scattering matrix. In 
this Good- Walker type formalism, diffractively produced 
hadrons at a given vertex are considered as a single 
hadronic state described by the wave function 'ilo, which 
is orthonormal to the wave function ^'/j of the incoming 
hadron, < >= 0. We introduce two wave func- 

tions 5*1 and ^'2 which diagonalize the 2x2 interaction 
matrix T 

Kk' *^*fc|T|*.- ^k' >= A,^kS^,^> Sk,k'- (II.l) 
In this representation the observed states are written 

*h = a*i+/3*2, (11.2) 

*I5 = +a*2, (11.3) 

where, + (3"^ ~ 1. 

Using Eq. (flLT|) we can rewrite the unitarity equations 

/mv4,,fc(s,6) = |A,,fc(s,5)|2 + G^:)_,(s,5), (II.4) 

where GY\, is the summed probability for all non diffrac- 
tive inelastic processes induced by the initial (i, k) states. 
The simple solution to Eq. pi.4[) has the form obtained 
in a single channel formalism 3], 



Ai^kis,b) = q 1 - exp - 



(H.5) 
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Gl^{s,b) = l-cxp{-n,Ms,b)). (II.6) 

From Eq. (jll.6[) we deduce the probability that the ini- 
tial projectiles (i, k) reach the final state interaction un- 
changed, regardless of the initial state re-scatterings, is 
given by Pfj, = exp (-^li,fc(s, b)). 

In general, we have to consider four possible (i, k) re- 
scattering options. For initial p-p (or p-p) the two quasi- 
elastic amplitudes are equal Ai^2 = ^2,1, and we have 
three re-scattering amplitudes. The corresponding elas- 
tic, SD and DD amplitudes are 

aeiis, b) = i{a^Ai^i + 2a2/32Ai,2 + /3^^2,2}, (11.7) 

asd{s, b) = ial3{-a^Ai^i + (a^ - (3'^)Ai,2 + P^A2,2}, 

(II.8) 

add = la^PHAi^i - 2Ai,2 + ^2,2}- (11.9) 

Adjusted parameters are introduced to obtain explicit 
expressions for the opacities ili,k{s,b). 

In the following we shall consider Regge and non Regge 
options for the dynamics of interest. We use a simple 
general form for the input opacities, 

n,,k{s,b) = i^,^Us)Tis,b). (11.10) 

'^^,k{s)=al, (^^y. (11.11) 

The input b-profiles F.^ ^ (s, b) are assumed to be Gaus- 
sians in b, corresponding to exponential differential cross 
sections in t-space, 

^^,kis,b) = J exp U-J^] , (n.12) 

RlAs) = Rh.k + '^Cln{s/so). (n.i3) 

Rq.i 2 = ^Rq i 1 and i?Q.2 2 = 0. Our parametrization 
is compatible with, but not exclusive to, a Regge type 
input. 

III. FITS AND PREDICTIONS 

We have studied three models, with different param- 
eterizations of ^li^k, which were adjusted to the ISR- 
Tevatron experimental data base, specified above. Note 
that the fit has, in addition to the contribution in the 
form of Eq. (|II.10|) . also a secondary Regge sector (see 
Ref. d, 01 ) • This is necessary, as the data base contains 
a relatively small number of experimental high energy 
measured values, which are independent of the Regge 
contribution. We do not quote the values of the Regge 
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TABLE I: Fitted parameters for Models A, B(l) and B(2). 



parameters, as the goal of this paper is to obtain pre- 
dictions in the LHC and Cosmic Rays energy range. At 
W=1800GeT^ the Regge sector contribution is less than 
1%. However, it is essential at the ISR energies. 

Model A is a simplified two amplitude version of the 
two channel model, in which we assume that add is small 
enough to be neglected. As such, this model breaks 
Regge factorization. The model was presented and dis- 
cussed in Ref.[3|. The parameters of Model A were ob- 
tained from a fit to a 55 experimental data points base 
and are listed in Table 1 with a corresponding /(d-o.f) 
of 1.50. Note that in Model A the (1,1) amplitude cor- 
responds to while the (1,2) amplitude corresponds 
to An = ni^i - ni^2- See Ref. [4]. 

Model B denotes our three amplitude model where the 
5 published DD cross section points [l3| are contained 
in the fitted data base. The three opacities are taken 
to be Gaussians in 6. If we assume the soft Pomeron 
to be a simple J pole, its coupling factorization implies 

= ^J^tl^<~^2■ We denote this Model B(l). The 

fit obtained is not satisfactory, with a x^/((i.o./.)=2.30. 

We have, also, studied Model B(2) in which coupling 
factorization is not assumed. Accordingly, a^i, ai2 and 
a2 2 are independent fitted parameters of the model. The 
model with a x^/ {d.o.f.) = 1.25, provides a very good re- 
production of our data base. In Model B(2) the leading 
t channel exchange is not a simple J pole. It is com- 
patible with a model[l2| we have suggested a while ago 
in which the soft Pomeron dominated photo and low 
DIS, is perceived as the saturated soft (low Q"^) Hmit 
of the hard Pomeron dominated (high Q^) hard DIS. A 
major deficiency of Model B(2) is that it predicts dips 
in at small t values, which are not observed experi- 
mentally. This problem is common to all eikonal models 
which assume Gaussian b-profiles. Consequently, Model 
B(2) is valid only in the narrow forward t cone, where it 
reproduces approximately 85% of the overall data very 
well. We shall discuss this problem in some detail in the 
Discussion Section. 

Model B(2) cross section and slope predictions at ultra 
high energies are summarized in Table 2. Note that R^i = 
(^eil'^tot and Rd = {(Jei + crdiff)/crtof At LHC (W=14 
TeV) our predicted cross sections are: atot — 110.5 m&. 
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TABLE 11: Cross sections and elastic slope in Model B(2). 



Uei — 25.3 m6, Usd — 11.6 m6 and add = 4:.9 mb. These 
predictions are slightly higher than those obtained^ in 
Model A. The corresponding forward slopes are: B^i = 
20.5Gey-2, B,d = 15.9 GeV-^ and Bdd = 13.5 GeV^-^. 
We calculate, also, p = 0.125. The calculations of Bg^, 
Bj^ii and p were executed with the fitted parameters of the 
model. For the record we have checked that we reproduce 
also the UA4, CDF and E710 Bgd and p data points. 



IV. SURVIVAL PROBABILITIES 

In the following we shall limit our discussion to the 
survival probability of Higgs production in an exclusive 
central diffractive process, calculated in our model. For 
a general review see Ref . |1| . 

In our model we assume an input Gaussian b- 
dependence also for the hard diffractive amplitude of in- 
terest. Its input, when convoluted with the soft (i,k) 
channel, is 



f^5 = ^5(^)rf,(&), 



,H 



(IV. 14) 



(IV. 15) 



(IV.16) 



The structure of the survival probability expression is 
shown in Fig. [TJ The corresponding general formulae for 
the calculation of the survival probability for diffractive 
Higgs boson production have been discussed in Refs.|[l|, 
[1,1^. Accordingly, 



(IV. 17) 
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FIG. 1: Survival probability for exclusive central diffractive 
production of the Higgs boson 
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(l-A5(.,(bi+b2)))}2, (IV18) 

D{s)^ J hi b2{ Ah {s,bi) Ah {s,b2)}^(W. 19) 

As denotes the soft strong interaction amplitude given by 
Eq. ( jlLll). Using Eg. (jUll- Eq. (HLQ]) . the integrands of 
Eq. PV.ISP and Eq. PV.IOP are reduced by eliminating 
common s-dependent expressions. 

Nis) = [ cPbid%2{{l-aei{s,b))APPibi)APPib2) 



-asd{s, b) (<(6i)A?r(62) + A«'(6i)<(62); 
-addis, b)A%\bM^^{b2)}\ (IV.20) 

j d%A{A^S(b,)A^^{b2)Y- (IV.21) 
Following Refs.P, Q we introduce two hard 6-profiles 
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H 



(IV.22) 
(IV.23) 



The hard radii Rf^. and cross section coefficients Vp-^p 
and Vp^d are constants derived from HERA J/^ elastic 
and inelastic photo and DIS production[l^. [ist (see, also, 
Ref.i). B^i = 3.6Gey-2, Bg = IGeV-^, Vp^p = VS 
and Vp^d = 1- have been taken from the experimental 
HERA data on J/^- production in HERA[i3,1i1]. 

Using Eq. (|IV.17|) -Eq. (|IV.21|) we calculate the survival 
probability for exclusive Higgs production in central 
diffraction. Sjj has been calculated[lj in the two ampli- 
tude Model A. The resulting Sj^ = 0.027 is essentially 
the same as the predictions of KMR[i3|. Our present re- 
sults, obtained in the three amplitude B Models, indicate 
a reduction of the output value of Sj^ . Its LHC value in 
Model B(l) is 0.02, and in Model B(2) it is 0.007. We 
note that, our Model B(l) result is compatible with the 
result of Ref.[l^. We shall return to this issue in the 
Discussion Section. 



4 



AMPLITUDE ANALYSIS 



The basic amplitudes of the GLM two channel model 
are Ai^i, Ai^2 and ^2,2, whose h structure is specified in 
Eq. (|II.5|1 ). These are the building blocks with which we 
construct a^i, agd and add (Eq. (jII.7P -Eq. (jll.9|l ). The 
Ai^k amplitudes are bounded by the black disc unitarity 
bound of unity. Checking Table 1, it is evident that in 
both Model B(l) and B(2) ^l2,2 is much larger than the 
other two fitted opacities. As a consequence, the ampli- 
tude ^2,2(5,^) reaches the unitarity bound of unity at 
low energies. Similarly, the output amplitude Ai_2(s,6) 
of Model A reaches unity at approximately LHC en- 
ergy. The observation that one, or even two, of our 
Ai^k{s,b)=l does not imply that the elastic scattering 
amplitude has reached the unitarity bound at these (s, b) 
values. aei{s,b) reaches the black disc bound when, and 
only when, Ai^i{s,b)=Ai^2{s,b)=A2^2{s,b)=l. In such 
a case we also obtain, that asd{s,b)=addis,b)=0. This 
result is independent of the fitted value of p. 

Model B(2) predictions of Uei over a wide range of en- 
ergies are presented in Fig. [2l A fundamental feature of 
Models A, B(l) and B(2) is that a^i approaches the black 
disc bound at 6 = very slowly, reaching the bound at 
energies higher than the GZK knee cutoff. If correct, 
this feature implies that aei does not reach the black disc 
bound over the entire accessible spectrum of Cosmic Rays 
energies, even though it gets monotonically darker. 

The explanation of this behavior, in our presentation, 
is simple. Checking the values of /3 and ^. correspond- 
ing to the 3 models (see Table 1), we note that fti^i is 
smaller by 1-3 orders of magnitude relative to 02,2 (^^1,2 
in Model A). The consequent ad can reach the black 
disc bound only when fii^i is large enough so that Ai^i 
approaches unity, fli^i grows slowly like W'^'^ (modulu 
InW). Hence, the slow approach of Uei toward the black 
disc bound. This result is incompatible with the output 
of Ref . in which Uei reaches the black disc bound ap- 
proximately at the LHC. In our presentation it implies 
that unlike our models, in the KMR model there is rela- 
tively small variance in the weights of the 3 components 
of the proton wave function. 

A consequence of the input fli^k being large at small 
6, is that Pfi^{s,b) is very small at 6 = and monoton- 
ically approaches its limiting value of 1, in the high b 
limit. As a result, given a diffractive (non screened) in- 
put, its output (screened) amplitude is peripheral in b. 
This is a general feature, common to all eikonal models 
regardless of their b-profiles details. The same is, true, 
also, with regard to diffractive Good- Walker channels, 
which are contained in f2i_/c. This implies a non trivial t 
dependence of dadif f{M^-^j)/dt in the diflfractive chan- 
nels. These qualitative features are induced by Model 
A, B(l) and B(2), even though their detailed behavior 
are not identical. Given the deficiencies of our b-profiles. 



500 TeV 




1 1.5 2 

b in fm 

FIG. 2: b dependence of a^i in Model B(2) at different energies 



we refrain from giving any specific predictions beside the 
general observation stated above. 

The general behavior indicated above becomes more 
extreme at ultra high energies, when aei continues to ex- 
pand and gets darker. Consequently, the inelastic diffrac- 
tive channels becomes more and more peripheral and rel- 
atively smaller when compared with the elastic channel. 
At the extreme, when a^^s^b) — 1, a^d — o-dd — 0. 
We demonstrate this feature and its consequence at the 
Planck mass in Fig. [21 As the black core of Uei expands, 
the difference between Models A, B(l) and B(2), con- 
sidered in this paper, diminishes, being confined to the 
narrow b tail where aei{s,b) < 1. The above observa- 
tions may be of interest in the analysis of Cosmic Ray 
experiments. 



VI. DISCUSSION 

It is interesting to compare our model and its out- 
put with a different eikonal model recently proposed by 
KMR 14 1 extending earlier versions ^iJl . The two mod- 
els were constructed with very similar objectives but are 
fundamentally different in their conceptual theoretical in- 
put, data analysis and output results. 
1) The input of KMR is a conventional Regge model in 
which high mass diffraction, initiated by Pomeron en- 
hanced diagrams, is included. GLM is a phenomenologi- 
cal parametrization in which we assume diffraction to be 
strictly Good- Walker type, with no high mass diffraction 
distinction. We formulate our input in a general form 
consistent with Regge, but not exclusively so. Our sta- 
tistically preferred non factorizeable Model B(2) is com- 
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patible with a partonic interpretation which considers the 
soft "Pomeron" to be a low high density Hmit of the 
hard Pomeron llj. The GLM "Pomeron" is not a Regge 
simple J-pole, it does not include Pomeron enhanced di- 
agrams, which are essential in the construction of KMR. 

2) Since multi-Pomeron vertices are included in KMR, 
they had to fix a' — 0. In order to maintain the exper- 
imentally observed forward f-cone shrinkage, they con- 
structed a high absorption eikonal model in which the 
input is non conventional A = 0.55. With this input, 
KMR obtain an approximate DL behavior in the ISR- 
Tevatron range. However, at higher energies their effec- 
tive A becomes monotonically smaller (its value in the 
Tevatron-LHC range is reduced to 0.04) which results 
in a very slow rise of a tot and (TeZ- GLM is a weak 
screening eikonal model. Its fitted input is A = 0.15 
and C = a' = 0.17. With this input, GLM total cross 
sections are compatible with DL over the wide ISR-GZK 
range. 

3) The goal of both GLM and KMR is to adjust the 
model parameters of their vacuum t exchange " Pomeron" 
input, so as to predict and calculate observables and fac- 
tors of interest at the LHC and Cosmic Rays. Both mod- 
els adjust more than 10 free parameters. Only CERN- 
UA4 and Tevatron energies are sufficiently high to justify 
neglecting the contribution of the secondary Regge sec- 
tor. This limited data base is not sufficient to adjust the 
"Pomeron" free parameters. GLM chose, therefore, to 
construct a model containing also the secondary Regge 
sector and fit the extended data base spanning the ISR- 
Tevatron energy range. KMR constrain their parameter 
adjustment to the small data base of the highest ener- 
gies. In our opinion the KMR procedure is not adequate. 
Indeed, their reconstruction of at the 3 highest avail- 
able energies is remarkably similar to a fit they made a 
few years ago with different parameters, notably a con- 
ventional A input. 

4) GLM and KMR determine their input opacities in 
completely different procedures which define their (dif- 
ferent) data bases. GLM approach is that a model which 
takes into account diffractive re-scatterings of the ini- 
tial projectiles has to reconstruct properly the diffractive 
cross sections, which are, thus, included in its fitted data 
base. KMR goal is to reconstruct aei{s,b) for which the 
diffractive components are needed. To this end they fit 

neglecting an explicit fit of the diffractive channels. 
Obviously, combining both GLM and KMR data bases is 
advisable. Regretfully, we were unable to obtain good si- 
multaneous reproduction of such an extended data base. 
The question, is thus, which model provides a better ap- 
proximation for the input opacities. 

5) The b-distributions of aei{s,b) in GLM are signifi- 
cantly different from KMR. GLM obtain a relatively wide 
b distribution compared with a narrower one in KMR. 
aei{s, 6 = 0) in KMR is consistently larger than in GLM, 
approaching the black disc bound much faster than in 



GLM. Regardless of these differences, the corresponding 
values of atot and aei in both models in the UA4- Tevatron 
range are compatible. Such compatibility can exist only 
over a relatively narrow energy band and it cannot per- 
sist over a wide energy range. Indeed, the two models 
have different LHC and Cosmic Rays predictions, which 
hopefully will be tested soon. Our inability to repro- 
duce outside the narrow forward t cone implies a 
deficiency in our aei at large b. We are not clear if this 
deficiency is reponsible for the small Sjj obtained in our 
Model B(2). Note, that even though our factorizable 
Model B(l) has the same feature of spurious dips outside 
the very forward at t cone, its predicted Sj^ is 0.02 which 
is compatible with KMR. 

6) In our opinion, the data adjustment procedure adopted 
by KMR are not adequate. Our approach is to quantify 
our fit by minimizing its x^. KMR reject any statisti- 
cal approach to their data analysis. They tune many 
of their parameters by eye and refrain from a quantified 
assessment of their output. The difference between the 
procedures adopted by the two groups is cardinal, as one 
is unable to make a systemic evaluation of the KMR out- 
put. 

7) The difference between the Sj^ predictions of GLM 
and KMR are intriguing and reflect the sensitivity of 
Sjj to each model input. Sj^ is calculated as a con- 
volution of the hard amplitude for Higgs production 
and the soft probability P/^,(s,6). The hard amplitude 
features needed for this calculation in our model are 
the hard slopes i?^, Bf^ and cross section coefficients 
Vp_^p Vp__^^, determined from the HERA measured ll.[l3l| 
in J/^ photo and DIS elastic and inelastic produc- 
tion. Our sensitivity to these parameters is shown in 
Fig. [3l Note that when we change the value of Bf^ , we 
keep the ratio V^^^/ B^ unchanged. Doing so we do 
not change the cross section of the reaction 7 -I- p — > 
J/* + X(M < 1.6 GeV). KMR calculation is simpler in 
as much as they consider just the elastic hard slope. In 
our opinion there is a gap between the sophistication of 
KMR soft model and the simplicity of their hard approx- 
imation. Since S'jj is obtained from a convolution of the 
two terms it is not clear what is the contribution of KMR 

hard term to the margin of error in their calculation of 
02 

KMR estimate their margin of uncertainty to be a 
factor of 2.5. Since our uncertainty derives from similar, 
though not identical, sources, our assessment is similar. 
As we saw, both GLM and KMR models are partially 
deficient. We noted that these are based on the different 
conceptual constructions and data analysis procedures 
of the two models. A discrimination between the two 
models depends on experimental results which are 
expected to become available within the next few years. 
In the following we list a few: 

1) GLM predictions for a tot and aei at the LHC are 
20% higher than the corresponding KMR values. This 



6 



4.5 

4 
3.5 

3 
2.5 

2 
1.5 

1 

0.5 






Bill 


\ S^(LHC)in% 


^ 0.5 GeV"^ 




, 1 GeV"^ 


■ .... 1 .... 1 .... 1 .... 1 ... . 


, 1.5 GeV"^ 



3.5 



4.5 



5 h5.5 , 6 

B"(GeV-) 



6.5 



FIG. 3: The dependence of at the LHC on Sff and B^„, 
the slopes for the hard cross sections. 



is a fundamental difference since the output energy 
dependence of GLM, which is a weak screening model, 
is compatible with an effective A = 0.08 all through the 
Tevatron-GZK energy range. In the KMR model the 
effective A is reduced rapidly due to the very strong 
screening which is inherent to this model. Hence, the 
KMR cross sections grow very moderately above the 
Tevatron energy. 

2) The difference between the two models becomes more 
distinguished at Cosmic Rays energies. This may be 
checked by the Auger experiments where we expect soon 
some cross section results at energies spanning up to W 
= 100-150 TeV. 

3) A basic feature particular to the KMR model is a 
contribution to diffraction which originates from the 
Pomeron induced diagrams which are not contained in 
GLM. As a result, both asd and add predicted by KMR 
are larger than GLM. These differences are very signif- 
icant for the DD channel where the KMR prediction 
at LHG is almost a factor of 3 larger than GLM. Note, 
that since diffraction in GLM is Good- Walker type, our 
predicted elastic and diffractive cross sections satisfy the 
Pumplin bound|17j|. aei{s,b) 



CFdiff{s,h) < ^atot- This 
bound does not aply to KMR, in which a significant part 
of its diffractive cross section originate from Pomeron 
enhanced contributions. 

4) An estimate of value can be obtained, at an early 
stage of LHC operation, through a measurement of the 
rate of central hard LRG di-jets production (a GJJG 



configuration) coupled to a study of its expected rate in 
a non screened pQCD calculation. 
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